Strong longitudinal color-field effects in 
pp collisions at energies available at the CERN Large Hadron Collider. 
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We study the effect of strong longitudinal color fields (SCF) in p+p reactions up to Large Hadron 
Collider energies in the framework of the HIJING/BB v2.0 model that combines (collinear factor- 
ized) pQCD multiple minijet production with soft longitudinal string excitation and hadronization. 
The default vacuum string tension, Ko = 1 GeV/fm, is replaced by an effective energy dependent 
string tension, «(s) = ko(s/so) ' 06 that increases monotonically with center-of-mass energy. The 
exponent A = 0.06 is found sufficient to reproduce well the energy dependence of multiparticle ob- 
servables in RHIC, Tevatron, as well as recent LHC data. This exponent is half of that predicted 
by the Color Glass Saturation (CGC) model, Xcgc = 0.115, where gluon fusion multiparticle pro- 
duction mechanisms are assumed. In HIJING/BB v2.0, the rapid growth of dN c h/dn with energy 
is due to the interplay of copious minijet production with increasing of strong color field (SCF) 
contributions. The large (strange) baryon-to-meson ratios measured at Tevatron energies are well 
described. A significant enhancement of these ratios is predicted up to the highest LHC energy (14 
TeV). The effect of J J loops and SCF on baryon-anti-baryon asymmetry, and its relation to baryon 
number transport, is also discussed. 

PACS numbers: 12.38.Mh, 24.85.+p, 25.40Vc, 25.75.-q 



I. INTRODUCTION 



With the commissioning of the Large Hadron Collider 
(LHC), it will soon be possible to test models of multipar- 
ticle production in hadron-hadron collisions up to an en- 
ergy of 14 TeV. Charged particle densities, dN c h/drj, es- 
pecially the values at mid-rapidity and their dependence 
on center-of-mass energy yfs (c.m.s.) are important for 
understanding the mechanism of hadron production and 
the interplay of soft and hard scattering contributions in 
the LHC energy range. The rate of parton-parton and 
multi parton-parton (MPI) scattering are strongly corre- 
lated to the observed particle multiplicity (related also to 
initial entropy and initial energy density generated in the 
collision process). New data on inclusive charged particle 
distributions from the LHC in pp collisions have become 
available (lHl^ . These results complement previous data 
on pp and pp collisions taken at lower energies yfs = 0.02 
- 1.96 TeV flli - Eijl ]. Many of these measurements have 
been used to constrain phenomcnological models of soft- 
hadronic interactions and to predict properties at higher 
energies [2~4H37l ]. The recent LHC data may lead to a 
better theoretical understanding based on a Quantum 
Chromodynamics (QCD) approach [38l - l4l| . 

The Heavy Ion Jet Interacting Generator (HIJING) 
[12 and HIJING/BB vl.10 models H3 have been used 
extensively to study particle production in pp collisions 
and to determine the physical properties of the ultra- 
dense matter produced in relativistic heavy-ion collisions. 
In the LUND Q and Dual-Parton (DPM) [H] models 
multi QCD strings or flux tubes were proposed to de- 
scribe soft multiparticle production in longitudinal color 
fields. Color exchange between high x partons in the 



projectile and target create confined color flux tubes of 
tension (sa 1 GeV/fm), that must neutralize through pair 
production or color singlet hadronization approximately 
uniformly in rapidity. In nucleus-nucleus collisions, the 
A 1 / 3 enhancement of the local parton density of high x 
partons allows for higher color Casimir representations to 
be excited. Those flux tubes with stronger longitudinal 
color fields than in average pp reactions have been called 
color ropes fdH ] and naturally have higher string tensions 
[13]. Recently, an extension of Color Glass Condensate 
(CGC) theory has prop osed a more detailed dynamical 
"GLASMA" model HI pi of color ropes. 

In the HIJING model (42| the soft beam jet fragmenta- 
tion is modeled by simple diquark-quark strings as in the 
LUND model with multiple gluon kinks induced by soft 
gluon radiations. Hard collisions arc included with stan- 
dard perturbative QCD (pQCD) as programmed in the 
PYTHIA generator [50|. However, HIJING differs from 
PYTHIA by the incusion of a geometric scaling multi- 
ple jet production model. Thus this model contains both 
longitudinal field induced soft beam jet multiparticle pro- 
duction and collinear factorized pQCD based hard mul- 
tiple jet production for pt > po — 2 GeV/c. 

A systematic comparison with data on pp and pp col- 
lisions in a wide energy range [42j revealed that minijet 
production and fragmentation as implemented in the HI- 
JING model provide a simultaneous and consistent expla- 
nation of several effects: the inclusive spectra at moder- 
ate transverse momentum (pt); the energy dependence of 
the central rapidity density; the two particle correlation 
function; and the degree of violation of Koba-Nielsen- 
Olesen (KNO) scaling [5l|, (HJ up to Tevatron energy 
(yfs = 1.8 TeV). However, the model failed to describe 
the dependence of the mean value of transverse momcn- 
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turn (< pt > ) o n charged particle multiplicity (iV c h). 
Wang argued (52] that by requiring high within a 
limited pseudorapidity (n) range one necessarily biases 
the data towards higher pt minijets, hence the observed 
increase of < pr > versus (vs.) dN c h/dr] (52J. This effect 
has also been associated with the presence of transverse 
flow of the hadronic matter (53l . |54| and was proposed 
as possibly due to quark-gluon plasma (QGP) formation 
already in pp collisions. 

Initial states of color gauge fields produced in high- 
energy heavy-ion collisions have also recently been dis- 
cussed in Ref. [5^]. Decay of a strong color electric field 
(SCF) (E > E cnUca i=10 18 V/m) due to the Schwinger 
mechanisms [56[ plays an important role at the initial 
stage of heavy-ion collisions at ultra-relativistic energies. 
A thermalization scenario based on the analogy between 
Schwinger mechanisms and the Hawking-Unruh effect 
has been proposed [57} • It was also suggested that the 
back-reaction and screening effects of quark and anti- 
quark pairs on external electric field could even lead to 
the phenomenon of plasma oscillations [58 60]. 

Recently, the Schwinger mechanism has been revisited 
[6~J ] and pair production in time-dependent electric fields 
has been studied Js3|. It was concluded that particles 
with large momentum were likely to have been created 
earlier, and for very short temporal widths (At 10t c , 
where the Compton time t c = l/m c ) and as a conse- 
quence the Schwinger formula could underestimate the 
reachable particle number density. In previous papers, 
we have shown that the dynamics of strangeness produc- 
tion in pp and Au + Au collisions at Relativistic Heavy 
Ion Collider (RHIC) energies deviates considerably from 
calculations based on Schwingcr-like estimates for homo- 
geneous and constant color fields and point to the con- 
tribution of fluctuations of transient strong color fields 

(scf) [i3H. 

In this paper we extend our study of the dynamic con- 
sequences of SCF in the framework of the HIJING/BB 
v2.0 model 64] to particle production in hadron-hadron 
collisions at LHC energies. We explore dynamical effects 
associated with long range coherent fields (i.e, strong lon- 
gitudinal color fields, SCF), including baryon junctions 
and loops [63l Wi\ . with emphasis on the new observ- 
ables measured in pp collisions at the LHC by the AL- 
ICE, ATLAS and CMS collaborations. Our study aims 
to investigate a broad set of observables sensitive to the 
dynamics of the collisions, covering both longitudinal and 
transverse degree of freedom. In addition, this study is 
intended to provide the pp baseline to future extrapo- 
lations to LHC studies of proton-nucleus (p + A) and 
nucleus-nucleus (A + A) collisions. 



mented. In a dual superconductor model of color con- 
finement for the three-quark positioning in a Y - geom- 
etry, the flux tubes converge first toward the center of 
the triangle and there is also another component that 
runs in the opposite direction (like a Mercedes star) [68j . 
Unlike the conventional diquark-quark implemented in 
LUND and the HIJING model jH, the HIJING/BB vl.10 
[43| model allows the diquark-quark to split with the 
three independent flux lines tied together with an eyfe 
a junction and terminating on three delocalized funda- 
mental Casimir quarks. We introduced [67j a new ver- 
sion (v2.0) of HIJING/BB that differs from HIJING/BB 
vf.IO [43[ in its implementation of additional more com- 
plex flux topologies via junction anti-junction ( JJ) loops. 
We parametrize the probability that a junction loop oc- 
curs in the string. Moreover, we enhance the intrinsic 
(anti)diquark-quark px kick (by a factor / = 3) of all 
(q-qq) strings that contain one or multiple JJ loops. The 
reason for this is the mechanism behind the dynamic of 
diquark-quark breaking (see Ref. [13] for details). 

In string fragmentation phenomenology, it has been 
proposed that the observed strong enhancement of 
strange particle production in nuclear collisions could be 
naturally explained via strong color field (SCF) effects 
[I?} . For a uniform chromoelectric flux tube with field 
(E), the pair production rate HE 53, EH) P er um t volume 
for a (light)heavy quark (Q) is given by 
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where Q 
(bottom) . 



qq (diquark), s (strange), c (charm) or b 



The current quark masses are m q 



0.45 



GeV |69(|, m s = 0.12 GeV, m c = 1.27 GeV, and m b = 
4.16 GeV pro( | . The constituent quark masses of light 
non-strange quarks are M,., a — 0.23 GeV, of the strange 
quark is M s =0.35 GeV (71|, and of the diquark is M qq = 
0.55 ±0.05 GeV 

Note that k = \eE\ e ff = ^C 2 {A)/C 2 (F) ko is the 
effective string tension in low energy pp reactions with 
Ko ~ 1 GeV/fm and Ci(A), C^F) are the second order 
Casimir operators (see Ref. (47[). An enhanced rate for 
spontaneous pair production is naturally associated with 
"strong chromo- electric fields", such that k/itVq > 1 at 
least some of the time. In a strong longitudinal color 
electric field, the heavier flavor suppression factor 'Jqq 
varies with string tension via the well known Schwinger 
formula |56J|, since 



1 qq 



exp 



m 2 q ) 



< 1 



Ko 



(2) 



II. OUTLINE OF HIJING/BB V2.0 MODEL. 

In HIJING/BB v2.0, in addition to conventional quark- 
diquark longitudinal electric fields, novel color flux topol- 
ogy junction anti-junction (JJ) loops are also implc- 



for Q = qq, s, c or b and q — u, d. In our calculations, 
we assume that Mjf = 0.5 GeV, Mf = 0.28 GeV, Mf 
= 1.30 GeV. Therefore, the above formula implies a sup- 
pression of heavier quark production according to u : d 
: qq : s : c w 1 : 1 : 0.02 : 0.3 : lO" 11 for the vac- 
uum string tension kq — 1 GeV/fm. For a color rope, 
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on the other hand, if the average string tension value k 
increases, the suppression factors Jqq increase. We show 
below that this dynamical mechanism improves consider- 
ably the description of the strange meson/hyperon data 
at the Tevatron and at LHC energies. 

Saturation physics is based on the observation that 
small-x hadronic and nuclear wave functions, and, thus 
the scattering cross sections as well, are described by the 
same internal momentum scale known as the saturation 
scale (Qsat)- A recent analysis of pp data up to LHC 7 
TeV has shown that, with the k? factorized (GLR) gluon 
fusion approximation [72| , the growth of the dN c h / drj can 
be accounted for if the saturation scale grows with c.m.s. 
energy as 



?o(s/*o) AcGC , 



(3) 



with Acgc ~ 0.115. The saturation scale is also increas- 
ing with atomic number as A 1 / 6 [5?]]. It was argued that 
the effective string tension («) of color ropes should scale 

withQL MM- 

However, in HIJING the string/rope fragmentation 
is the only soft source of multiparticle production and 
multiple minijets provide a semi-hard additional source 
that is computable within collinear factorized standard 
pQCD with initial and final radiation (DGLAP evolu- 
tion (zH). In order to achieve a quantitative descrip- 
tion, within our HIJING/BB framework we will show 
that combined effects of hard and soft sources of multi- 
particle production can reproduce the available data in 
the range 0.02 < yfs < 20 TeV only with a reduced de- 
pendence of the effective string tension on s/s. We find 
that the data can be well reproduced taking 

K{s) = K (s/ So ) a06 GeV/fm » Q Q Q sat (s), (4) 

where kq = 1 GeV/fm is the vacuum string tension value, 
so = 1 GeV 2 is a scale factor, and Qo is adjusted to give 
K = 1.88 GeV/fm at the RHIC energy y/s = 0.2 TeV. Our 
phenomenological k(s) is compared to Qs at (s) in Fig. [IJ 
where k = 1.40 GeV/fm at y/s = 0.017 TeV increases to 
k = 3.14 GeV/fm at Vs = 14 TeV. 

The energy dependence of the string tension leads to 
a variation of the diquark/quark suppression factors, as 
well as the enhanced intrinsic transverse momentum kx- 
These include i) the ratio of production rates of diquark- 
quark to quark pairs (diquark-quark suppression factor), 
7 qq = P(qqqq)/P(q , g); ii) the ratio of production rates 
of strange to non-strange quark pairs (strangeness sup- 
pression factor), 7 S = P(ss)/P(qq); hi) the extra sup- 
pression associated with a diquark containing a strange 
quark compared to the normal suppression of strange 
quark (j s ), j us = (P(usus)/P(udud))/(7 S ); iv) the sup- 
pression of spin 1 diquarks relative to spin ones (apart 
from the factor of 3 enhancement of the former based 
on counting the number of spin states), 710; and v) 
the (anti)quark (o'^ = \J k/kq ■ a q ) and (anti)diquark 

("q'q = \/ K / K o ' / ' °qq) Gaussian width. As an example 
we plot in Fig. [5] the energy dependence of the suppres- 
sion factor 7 S = s/u, when the string tension values k(s) 
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FIG. 1: (Color online) Energy dependence of the effective 
string tension k(s) in pp collisions (Eq. [4j . The values Qsa t 
(Eq.[3j) are from CGC model fits to LHC data from Ref. 0. 



are taken from Eq. 0] in comparison with the values pre- 
dicted using Q* at (Eq.[3]from the CGC model fit (2j|). 
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FIG. 2: (Color online) Energy dependence of the strange to 
light quark suppression factors, 7s = s/u, using k(s) from 
Eq. [4] and using Qs at (s) from Ref. [2!| are compared. (Note 
that this Schwinger suppression factor is not used in the CGC 
model). 

The contributions of multiple jets to the multiplicity 
distributions in pp and pp collisions have been studied in 
detail in Ref. [74[. Within the HIJING model, one as- 
sumes that nucleon-nucleon collisions at high energy can 
be divided into soft and hard processes with at least one 
pair of jet production with px > Po- A cut-off scale po 
in the transverse momentum of the final jet production 
has to be introduced below which the interaction is con- 
sidered non-pcrturbative and can be characterized by a 
finite soft parton cross section <7 so ft. The inclusive jet 
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cross section Oj Ct at leading order (LO) [75| is 



(Tjet 



s/i 



dp^dyidy 2 



1 da 



jet 



IT, 



2 dp|dyidj/2 ' 



(5) 



where, 



da 



jet 



da ab (s,t, u) 
dl 



(6) 

depends on the parton-parton cross section a ab and par- 
ton distribution functions (PDF) fa{ x iPr)- The summa- 
tion runs over all parton species; y\ and yi are the ra- 
pidities of the scattered partons; x\ and xi are the light- 
cone momentum fractions carried by the initial partons. 
The factor K w 2 accounts for the next-to-leading order 
(NLO) corrections to the leading order jet cross section. 
In the default HIJING model (v. 1.383), the Duke-Owens 
parameterization [76| of PDFs in nucleons is used. With 
the Duke-Owens parameterization of PDFs, an energy 
independent cut-off scale po =2 GeV/c and a constant 
soft parton cross section a so ft = 57 mb are sufficient to 
reproduce the experimental data on total and inelastic 
cross sections and the hadron central rapidity density in 
p+p(p) collisions [UHIl- Our results have been obtained 
using the same set of parameters for hard scatterings as 
in the latest version of HIJING (vl.383). 



III. CHARGED PARTICLES 
A. Charged Hadron Pseudorapidity 

Charged hadron multiplicity measurements are the 
first results of the LHC physics program. Data for pp col- 
lisions were reported by the ALICE, ATLAS, and CMS 
collaborations [H-tlH- The new data at mid-rapidity for 
non single diffractive interactions (NSD) and inelastic 
scattering (INEL) are shown in Fig. [3J which includes 
also similar results at lower energies. The main result is 
an observed sizeable increase of the central pseudorapid- 
ity density with c.m.s. energy 

The main contribution to the multiplicity comes from 
soft interactions with only a small component originat- 
ing from hard scattering of the partonic constituents of 
the proton. In contrast to the higher pt regime, well de- 
scribed by pQCD, particle production in soft collisions 
is generally modeled phcnomcnologically to describe the 
different pp scattering processes: elastic scattering (el), 
single diffractive (SD), double diffractive dissociation 
(DD), and inelastic non-diffractive scattering (ND). Ex- 
perimentally, minimum bias events are a close approxi- 
mation of NSD interactions, i.e., ctnsd = ftot — °ci — o"sd , 
where Ctot is the total cross section. The selection of 
NSD events is energy dependent and differs somewhat 
for different experimental triggers. The event selection 
of inelastic processes (INEL) includes SD interactions: 
cinel = otot — cr e i. The data therefore must be corrected 



for the SD component, involving model dependent calcu- 
lations. 

The results of the model calculations for both NSD 
(left panel) and INEL (right panel) are also shown in 
Fig. [3J Solid lines depict results including the strong 
color field (SCF) effects whereas the results without SCF 
effects are shown as dashed lines. Without SCF effects 
the model strongly overestimates the central charged par- 
ticle density. The absolute value and energy increase of 
the central rapidity density are well reproduced assuming 
the energy dependent string tension given in Eq. |4j At 
higher LHC energies (2.36 and 7 TeV), a discrepancy of 
10-15 % is observed. 
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FIG. 3: (Color online) Comparison of HIJING/BB v2.0 pre- 
dictions for central charged particle pseudorapidity density in 
pp and pp interactions for non-single-diffractive (NSD) (left 
panel) and inelastic (INEL) (right panel) interactions as a 
function of c.m.s. energy. The solid and dashed lines are 
the results with and without SCF, respectively. The data 
are from Refs. flHl, !, OH LIML1 [2l| (left panel) and from 
Refs. [3l. Iigl I2H |2^| (right panel). Only statistical error bars 
are shown. The open stars at 7 and 14 TeV (ALICE fit), are 
obtained by a power law fit to lower energy data from Refs. 

1, M. 



As the colliding energy increases, the rate of multi- 
ple parton interactions (MPI) also increases, producing a 
rise in the central multiplicity. The increase with energy 
in our phenomenology is due to the interplay of the in- 
creased mini-jet production in high colliding energy with 
SCF effects. For an increase of strangeness suppression 
factors due to an increase of string tension with energy 
(k = k(s)), the model predicts a decrease of produced 
pions due to energy conservation. Lower values of n(s) 
imply smaller values for strangeness suppression factors, 
therefore a higher multiplicity of mesons (mostly pions). 

Changing the effective value k(s) = 2.89 GeV/fm to 
k(s) = 2.0 GeV/fm results in an increased multiplicity of 
11% at 7 TeV where the effect is greatest. In addition, 
the multiplicity depends also on the value of the cut- 
off parameter pq. Low values of po imply high rates of 
parton-parton scattering and hence high levels of particle 
multiplicity. Evidently for increasing values of po the 
opposite is expected. At 7 TeV, where the effect is also 
the largest, changing po by ± 0.5 GeV results in a change 
of only =F 1.8% for the central pseudorapidity density. 
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FIG. 4: (Color online) Comparison of HIJING/BB v2.0 pre- 
dictions for the pseudorapidity distribution of charged parti- 
cles in p + p (p) collisions at various c.m.s. energies. The 
solid histograms are the results with SCF and JJ loops. The 
data are from Refs. [H (UA5), [TJ] (STAR), [3 (CDF), 
HH (ALICE), 0,0] (CMS). Only statistical error bars are 
shown. 



Data on the charged particle pseudorapidity distribu- 
tion are also available over a limited r] range 0[2| [!, H, [|| . 
These are presented in Fig. [4] where we also include re- 
sults obtained at 0.2 TeV by the UA5 [gland STAR 
fl3j collaborations, and with CDF results [l7| obtained 
at the Tevatron for pp collisions at 1.8 TeV. Consistent 
with the discussion above, a scenario with SCF effects 
(solid histograms) reproduces the measured multiplicity 
distributions well. At all energies considered, theoretical 
calculations predict a central dip at mid-rapidity that is 
consistent with the observations. At 0.9 and 2.36 TeV, 
the shape of the distribution measured by the ALICE 
collaboration is very well reproduced, while the CMS re- 
sults show a much flatter distributions than calculated 
(also the case at 7 TeV). Data over a larger rapidity 
range are needed to determine the shape of the falling 
density in the fragmentation region. For completeness, 
predictions at 10 and 14 TeV, the higher LHC energies, 
are also shown. 



B. Transverse momentum spectra 

The measured transverse momentum distributions for 
NSD events over an energy range y/s = 0.63 — 7 TeV 
are shown in Fig. [SJ These recent measurements are per- 



formed in the central rapidity region and cover a wide px 
range (0.15 < px < 10 GeV/c), where both hard and soft 
processes are expected to contribute. The data of ATLAS 
and CMS are measured in larger pseudorapidity intervals 
(|?7| < 2.5). In contrast, ALICE and CDF measurements 
are in a very central region (|?y| < 0.8 and \ri\ < 1.0, re- 
spectively) . The calculation takes into account the differ- 
ence in acceptance, but as can be concluded from Fig. 21 
this difference in pseudorapidity range has a negligible 
effect on the measured cross section. 



> 

u 





p- 

H 10 

5. 10 
5? 10 
S 10 

i 10 



10 
10 
10 
10 
10 



+ p, BB v2.0, charged, NSD 

■ 14TeV(X 100) 

■ 7 TeV (X 1 00) 

■ 2.36 TeV(X 10) 

■ 0.9 TeV 

0.63 TeV (X 0.01) 



*CMS, 7 TeV 
iCMS, 2.36 TeV 
■ CMS, 0.9 TeV 
t ATLAS, 0.9 TeV 

* ALICE, 0.9 TeV 

• CDF, 0.63 TeV 




****«. 

k = k(s) GeV/fm # *^ 



8 10 12 

p T [GeV/c] 



14 



FIG. 5: (Color online) Comparison with data of HIJING/BB 
v2.0 predictions of charged-hadron transverse momentum dis- 
tributions at LHC energies. The calculated spectra include 
the combined effects of SCF and JJ loops. The histograms 
and data have been scaled for clarity by the factors indicated. 
The data are from Refs. [f3 (CDF), Q (ALICE), 0] (AT- 
LAS), 0,0] (CMS). Error bars include only the statistical 
uncertainties. 



The model calculations including SCF effects describe 
the data well at ^/s = 0.63 TeV, but lead to a harder 
spectrum that observed at higher energy. The model 
gives a fair description of the spectral shape at low px 
but overestimate the data at high p T . The discrepancy 
is highest (up to a factor of three) at y/s = 7 TeV. We 
do not understand the source of this discrepancy and 
await higher px data to draw a firm conclusion. How- 
ever, in our phenomenology this could indicate that jet 
quenching, i.e. suppression of high px particles like that 
observed at RHIC energies in nucleus-nucleus collisions, 
could also appear in pp collisions in events with large mul- 
tiplicity. This overcstimation of high px yield leads also 
to a similar overcstimation of the mean transverse mo- 
mentum (< pt >) and of the correlation of mean < px > 
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as a function of N c h, which we do not discuss here. 

Within our model we generate events with different 
numbers of mini-jets up to some maximum values. High 
numbers of mini-jets lead to higher multiplicity events. 
For events with ten mini-jets (the maximum assumed in 
our calculation), the central charged particle pseudora- 
pidity density could increase up to f=a 20 and the total 
multiplicity could be greater than 150 at Tevatron energy 
(1.8 TeV). The measurements of two particle correlations 
over the entire azimuth could reveal the jet structure re- 
lated to high pt particles 78j. The study of these cor- 



relations in events with high multiplicity could help us 
to draw a firm conclusion with regard to a possible jet 
quenching phenomenon in pp collisions. 



IV. IDENTIFIED PARTICLE SPECTRA AND 
RATIOS 

A. Baryon-to-Meson ratio 

The pp single particle inclusive px spectra measure- 
ments are important for understanding collision dynam- 
ics, since the various particles show different systematic 
behavior, as observed at RHIC energy [78j]. Detailed the- 
oretical predictions for single inclusive hadron produc- 
tion (including strange hyperons) are discussed in this 
section. Baryon-to- meson ratios (B/M) are experimental 
observables that can be used at the LHC for investigat- 
ing multi-parton interactions and helping to understand 
the underlying physics [l§-[8l|. 



(a) p + p , s'£ = 1 .80 TeV 



+ 1 
°< 




^ (b) p + p , s'£ = 1 .80 TeV 


- *J* 




\/ m CDF data 




: k = *(s) GeV/fm 


K. = 1 .0 GeV/fm 

' 



p T [GeV/c] 



p [GeV/c] 



FIG. 6: (Color online) HIJING/BB v2.0 predictions at ^/s 
= 1.8 TeV of baryon and mesons transverse momentum at 
mid-rapidity (— 1 < y < 1) (a) and their ratios for single 
strange particles (A + A ) /Kg (b) in minimum bias pp col- 
lisions. The solid and dashed lines have the same meaning 
as in Fig . [3] Experimental results at y/s = 1.8 TeV are from 
Ref. [86| (CDF Collaboration). Error bars include only the 
statistical uncertainties. The ratios (b) have been calculated 
by us, dividing the spectra in the panel (a). 

Unexpectedly high B/M ratios observed in A + A colli- 
sions have been discussed in terms of recombination and 
coalescence mechanisms fo^l - l&Ij ] . Such high ratios at in- 
termediate pt were also reported in pp collisions at RHIC 



|85j and at the Tevatron [86[. In pp collisions, however, 
a coalescence/hadronization scenario is not favored due 
to low phase space density in the final state. Our HI- 
JING/BB model, with SCF effects included, provides an 
alternative dynamical explanation of the heavy-ion data 
at RHIC energies. We have shown that the model also 
predicts an increasing yield of (multi) strange particles, 
thereby better describing the experimental data ffrij . 
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FIG. 7: (Color online) Predictions of the HIJING/BB v2.0 
model for baryon and mesons transverse momentum in the 
rapidity range —2 < y < 2 (a) and their ratio (b) for minimum 
bias pp collisions at ^fs — 7.0 TeV. The solid and dashed 
lines have the same meaning as in Fig. O The experimental 
results (a) are from Ref. [871 (CMS Collaboration) . Error bars 
include only the statistical uncertainties. The ratios (b) have 
been calculated by us, dividing the spectra in the panel (a) . 

Figure [6] (b) shows a comparison of model predictions 
with CDF experimental data [86[ of the strange baryon- 
to-meson ratio (A + A°)/K% at 1.8 TeV. The particle 
Pt spectra are shown in the panel (a). The measured 
ratio is fairly well described within our phenomenology. 
The larger string tension parameterization results in a 
predicted increase of the ratio (A° + A )/ Kg by a fac- 
tor of « 10 at the Tevatron energy. The pt spectra 
show that the increased ratio is due almost entirely to 
an increase of the A cross section that is well described 
with k(s) = kq (s/so) ' 06 GeV/fm. The model underes- 
timates the kaon production by 15-25 %. 

The model predictions at 7 TeV, currently the maxi- 
mum energy where there are data [87J , are shown in Fig. [7] 
(b) . The model gives a good description of hyperon pro- 
duction (A + A ), for which an increase by a factor of 
ten is still predicted if SCF effects are considered. How- 
ever, our calculations underestimate by approximately a 
factor of two the yield of K s and, as a consequence, re- 
sult in a higher ratio than that observed. This result 
needs further investigation (theoretical and experimen- 
tal) on heavy flavour production at this energy. We note 
that the models PYTHIA [H, H3 and Energy-conserving 
Partons Off-shell remnants and Spliting of partons lad- 
ders (EPOS) [9(1 HH cannot reproduce the observed high 
B/M ratio (see Fig. 6 and Fig. 7 from Ref. |oj]). The 
preliminary data reported by the CMS collaboration in- 
dicate high hyperon yield. Comparisons with PYTHIA 
results show that this model significantly underestimates 
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the hyperon yields in pp collisions at 0.9 and 7 TeV [871 ]. 

The strange particle ratios could also be the manifes- 
tation of new collective phenomena. In the EPOS model 
such an increase is obtained if the production of a mini- 
plasma is considered in pp collisions 24 1, [jnj. If con- 
firmed by future measurements, the study of these ob- 
scrvables could open a perspective on new physics in pp 
interactions. 

Similar conclusions are obtained from the study of the 
proton/pion (p/tt + , v/tt~) ratios where data exist at 
lower energies [92, |93[ . These data, shown in Figure [31 
are limited to low px < 2 GeV/c. Adding SCF effects 
results in a very sizable increase of the ratio and our cal- 
culations provide a good description of the data in the 
measured range. However, as the calculations indicate, 
to draw a final conclusion measurements at intermediate 
and high pt are needed. 
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FIG. 8: (Color online) Comparison of HIJING/BB v2.0 pre- 
dictions with data on the non-strange baryon over meson ra- 
tios from minimum bias events in the rapidity range \y\ < 2. 
The solid and dashed lines have the same meaning as in Fig. [3] 
Experimental results for \y\ < 2 at yfs — 0.54 TeV (left panel) 
and at ^/s = 1.8 TeV (right panel) are from Ref. [H] (E735 
Collaboration). The results at mid-rapidity at yfs = 0.54 TeV 
(left panel) are from Ref. (UA2 Collaboration). Error 
bars include only the statistical uncertainties. 

To the extent that the LHC experiments are able to 
identify hadron species, such data will provide vital input 
to validate this interpretation. The model predictions at 
LHC energies for the px dependence of the p/ir~ ratio are 
shown in Fig. |H1 An enhancement up to the highest LHC 
energy and a weak energy dependence, with a saturation 
that sets in for a c.m.s. energy y/s > 2.36 TeV, is pre- 
dicted. Note that preliminary data at 0.9 TeV reported 
by the ALICE collaboration for px spectra of pions (tt + ) 
and protons (p) (94| cover to pt < 2.5 GeV/c. The model 
results, with SCF effects included (dot-dashed histogram 
in Fig. [9]) arc consistent with a p/tt + ratio derived from 
the spectra reported by ALICE at 0.9 TeV in Ref. [H]. 

In our approach, the dynamical mechanism that leads 
to such high values of B/M ratios is SCF appear- 
ing at the initial stage of the interaction. The SCF 
mechanism strongly modifies the fragmentation processes 
(strangeness suppression factors) and thus results in a 
huge increase of (strange)baryons. This interpretation is 
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FIG. 9: (Color online) Predictions of the HIJING/BB v2.0 
model of non-strange baryon over meson ratios (p/n~) for 
minimum bias events at mid-rapidity at LHC energies. The 
upper curves correspond to calculations which include the ef- 
fects of SCF and JJ loops. The lower curves corresponds to 
calculations without SCF effects. The results are at 0.9 TeV 
(dot dashed-histograms), at 7 TeV (dashed histograms), at 10 
TeV (dotted histograms) and at 14 TeV ( solid histograms). 



also supported by more sophisticated theoretical calcu- 
lations, in a scenario in which a time-dependent pulse 
for the initial strength of the color field is considered 
[95I |96| . The large enhancement of the baryon-to- meson 
ratios demonstrates that SCF could play an important 
role in multiparticlc production in pp collisions at LHC 
energies and that high energy density fluctuations can 
reach very high densities, potentially comparable to those 
reached in central Au + Au collisions at RHIC energies 

M. 



B. Baryon-anti- baryon asymmetry 

From the study of the baryon-anti-baryon asymme- 
try one can learn about the mechanism of baryon num- 
ber transport. Baryon production via the conventional 
default quark-diquark mechanisms in the Lund string 
model are known to be inadequate even in e + + e~ phe- 
nomenology. This is one of the main reasons for our con- 
tinued exploration of alternative baryon junction mech- 
anisms. The details of the new implementation of JJ 
loops in HIJING/BB v2.0 are described in Ref. 
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II A. In HIJING/BB v2.0 the main two mechanisms for 
baryon production are quark-diquark (q — qq) string frag- 
mentation and junction anti-junction (JJ) loops [68( in 
which baryons are produced approximately in pairs. In 
a junction loop a color flux line splits at some intermedi- 
ate point into two flux lines at one junction and then the 
flux lines fuse back at an anti-junction somewhere fur- 
ther along the original flux line. The distance in rapidity 
between these points is chosen via a Rcgge distribution 
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[64|. We assume that, out of the non single diffractive 
nucleon- nucleoli (NN) interaction cross section (ctnsd)> a 
fraction fjj = (Tjj/ (oinel — 0sd) of the events excite a 
junction loop. The probability that the incident baryon 
has a JJ loop in p(A)+A collisions after rihits simulated 
binary collisions is given by 



jj 



1 - (1 - /. 



jj) 



(7) 



where Ojj = 17 mb, ctsd is parametrized in the model, 
and the total inelastic nucleon-nucleon cross sections are 
calculated These cross sections imply that a junc- 
tion loop occurs in pp collisions with a rather high prob- 
ability (at RHIC energy fjj w 0.5 for ct inel = 42 mb). 
Taking a constant value for CTjj results in a decrease with 
energy of the probability Pjj, due to a faster increase 
with energy of ctinel relative to ctsd ■ The actual proba- 
bility is modified also by string fragmentation processes 
for which we consider a threshold cutoff mass M c = 6 
GeV/c 2 in order to have enough kinematic phase space 
to produce BB pairs. We have investigated the sensitiv- 
ity of the results to the value of parameters CTjj and M c 
and found no significant variation on pseudorapidity dis- 
tributions of charged particles for 15 mb < ct jj < 25 mb 
and for 4 GeV/c 2 < M c < 6 GcV/c 2 . 

Baryon number transport is quantified in terms of the 
rapidity loss {6yi oss = 6y hcam - 6y hal y 0n , where y bc am and 
2/baryon are the rapidity of incoming beam and outgoing 
baryon, respectively) and has been discussed within our 
model phenomenology for A + A collisions in Refs. (63l — 
|65| . It was shown that HIJING/BB vl.O overestimate the 
stopping power and give a mild energy dependence of 
net-baryons at mid-rapidity. The energy dependence of 
net-baryons at mid-rapidity per participant pair within 
HIJING/BB v2.0 is proportional to (s/s )~ 1/4+A/2 , sim- 
ilar to the dependence predicted in Rcf. [98| . with the 
assumption that JJ is the dominant mechanism. This 
dependence is obtained assuming the following parame- 
ters [65j|: so = 1 GcV 2 (the usual parameter of Regge 
theory), a(0) = 1/2 (the reggcon intercept of the tra- 
jectory), and ap(0) = 1 + A (where A « 0.01) for the 
pomeron intercept. 

Recently, the ALICE Collaboration reported results 
[9!| on mid-rapidity anti-proton-to-proton ratio in pp 
collisions at yfs = 0.9 and 7 TeV and equivalcntly the 
proton-anti- proton asymmetry, A = (N p — N p )/(N p + 
N p ). These data could be used to constrain Regge in- 
spired model descriptions of baryon asymmetry. The au- 
thors state that, within statistical errors, the observed 
p/p ratio shows no dependence on transverse momentum 
or rapidity in the limited measured acceptance (—0.8 < 
y < 0.8; < p T < 1 GcV/c). In Fig. [ID] (a) are com- 
pared the HIJING/BB v2.0 model predictions with the 
published data at 0.9 TeV. Our model predicts negligible 
dependence on pr for < pr < 2 GeV/c, and a slight 
Pt dependence at higher pr, where both effects (JJ loops 
and SCF) could contribute. The full calculation is shown 
by the solid histogram. The dotted line is the prediction 
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FIG. 10: (Color online) HIJING/BB v2.0 predictions p T dis- 
tributions (a) at mid-rapidity for the p/p ratio at yfs = 0.9 
TeV; rapidity distributions are shown in part (b). The re- 
sults are from three possible scenarios: without contributions 
from SCF and J J loops (dotted histograms), including only 
the effect of JJ loops (dashed histograms) and including both 
effects (SCF and J J loops) (solid histograms). The data are 
from Ref. [9!| (ALICE Collaboration). Only statistical errors 
are shown. 



without JJ loops and SCF effects while the dashed line 
includes the effect of JJ loops only. Over the measured 
ranges the rapidity distribution of the ratio p/p distribu- 
tion is not sensitive to the various scenarios presented. 
The scenario with combined effects results in a wider ra- 
pidity distribution at y > 3 (b). The narrow structure 
observed near y = has no physical significance and we 
believe that it is likely due to a numerical artifact of our 
current implementation of fragmentation scheme. 

Separate proton and anti-proton rapidity distributions 
are, however, much more sensitive to SCF effects, as seen 
in Fig. [TT] The model predicts a substantial increase (by 
a factor of w 5) for p(p), when SCF are taken into ac- 
count. Due to the high cutoff mass M c = 6 GeV/c 2 the 
effect of Jj loops is very small over the entire rapidity 
region. Analysis of the 7 TeV data leads to similar con- 
clusions, but shows less effect at high px- 
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FIG. 11: (Color online) Predictions of the HIJING/BB v2.0 
model for rapidity distributions of protons (a) and anti- 
protons (b) at y/s — 0.9 TeV. The histograms have the same 
meaning as in Fig. [TUJ 

Over the measured range the fccd-down-corrected data 
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for the p/p ratio rises from 0.957 ± 0.006 at 0.9 TcV to 
0.991 ±0.005 at 7 TeV j9|. Although the measured mid- 
rapidity ratio is close to unity there is a small but signifi- 
cant excess of protons over anti-protons corresponding to 
an asymmetry of A = 0.022±0.003 and A = 0.005±0.003 
at y/s = 0.9 and 7 TeV, respectively. Within our model 
using a Rcggc intercept of a$ = 1/2 results in a proton- 
anti-proton asymmetry A modol = 0.032 (0.002) at 0.9 TeV 
(7 TeV). The result overestimates the value of asymme- 
try at 0.9 TeV. These values show that a fraction of 
the baryon number associated with the beam particles 
is transported over rapidity intervals of less than seven 
units. However, a better understanding of the baryon 
transport would require measuring the rapidity depen- 
dence of the asymmetry over a larger range. 

In Ref. [64[ we discuss the predicted rapidity correla- 
tion length for production of baryon-anti-baryon pairs 
within our model. The predicted rapidity correlation 
length (1 — (^(O)]" 1 depends upon the value of the Regge 
intercept q(0) [43||. A value of a(0) ~ 0.5 leads to rapid- 
ity correlat ions in the range \db — y§\ ~ 2, while a value 
a(Q) ~ 1.0 [HI is associated with infinite range rapidity 
correlations. This kind of analysis and future measure- 
ments at LHC energies will help us to determine bet- 
ter the specific parameters characterizing JJ loops used 
in the calculations of baryon number transport and/or 
baryon-anti-baryon asymmetry while helping to establish 
the validity of Regge inspired models. 

V. SUMMARY AND CONCLUSIONS 

We have studied the influence of strong longitudinal 
color fields and of possible multi-gluon dynamics (gluon 
junctions) in particle production in pp collisions, with 
a focus on RHIC, to Tevatron, and LHC energies. We 
have investigated a set of observables sensitive to the dy- 
namics of the collisions, covering both longitudinal and 
transverse degree of freedom. A detailed comparison 
with newly available experimental data from the LHC 
has been performed. 

We found that the inclusion of the multiple minijet 
source limits the growth of the string tension k(s) to 
be approximately only linear as a function of saturation 
scale (5 sa t (Eq. [f]) , in contrast to recent approaches (29j 
where k(s) scales as Q^ at (Eq. [3]). The reason for this is 
that in the CGC model the collinear factorized minijet 
mechanism is suppressed by geometric scaling to much 
higher pt- Future measurements at LHC energies (y/s = 
7 and 14 TeV), extended to high pr, will help to clarify 
the validity of this mechanism. 

We have shown that SCF could play an important role 
in particle production at mid-rapidity in pp collisions. 
Our calculations show that high-energy density fluctua- 
tions in pp collisions at LHC can reach densities compa- 
rable to those reached in central nuclear (A+A) collisions 
at RHIC. A large enhancement of the (strangc)baryon- 
to-meson ratios that persists up to the highest LHC 
energy can be explained as an effect of SCF that ap- 



pears at the initial stage of the interaction. The mecha- 
nisms of hadron production are very sensitive to the early 
phase of the collisions, when fluctuations of the color field 
strength are highest. Strong Color Field effects arc mod- 
eled by varying the effective string tension that controls 
the qq and qqqq pair creation rates and strangeness sup- 
pression factors. SCF, therefore, may modify the frag- 
mentation processes with a resultant huge increase of 
(strange)baryons. 

We show that both JJ loops and SCF effects could 
play an important role in baryon (anti-baryon) produc- 
tion at mid-rapidity in pp collisions at LHC energies. In- 
troducing a new JJ loop algorithm in the framework of 
HIJING/BB v2.0 leads to a consistent and significant im- 
provement in the description of the recent experimen- 
tal results for proton-anti-proton and for baryon-anti- 
baryon asymmetry in comparison to the older versions 
HIJING/B or HIJING/BB vl.O 0. We have shown 
that baryon number transport is suppressed for 5y > 7, 
a result that is confirmed by recent ALICE measurements 

The present study is limited to the effect of initial 
state baryon production via possible junction dynamics 
in strong fields. It would be very useful to consider a 
generalization of back reaction effects [58| to the case 
not only of pair production relevant for mesons but to the 
more difficult three string junction configurations needed 
to describe baryon production. 

A greater sensitivity to SCF effects is expected also for 
open charm and bottom production [66j. In particular, 
measurements of rapidity and pr distributions for parti- 
cles involving charm and bottom quark, would provide an 
important test of the relevance of SCF fluctuations, help- 
ing us to determine values of the suppression factors Jqq 
(where Q = qq, s, c, 6), which have strong dependence on 
the main parameters of QCD (the constituent and cur- 
rent quark masses) and on the system size. Even though 
the success of this procedure has been clearly illustrated 
here, a fuller understanding of particle production and es- 
pecially of (multi)strangc particles in ultra-rclativistic pp 
collisions at the LHC remains an exciting open question, 
and will continue to challenge many theoretical ideas. 



VI. ACKNOWLEDGMENTS 

Acknowledgments: We thank S. Das Gupta and S. 
Jeon for useful discussions and continued support. We 
thank P. Levai for helpful discussions and suggestions 
throughout this project. VTP acknowledges computer 
facilities at Columbia University, New York, where part 
of these calculations were performed. This work was sup- 
ported by the Natural Sciences and Engineering Research 
Council of Canada. This work was also supported by the 
Division of Nuclear Science, of the U. S. Department of 
Energy under Contract No. DE-AC03-76SF00098 and 
DE-FG02-93ER-40764. 



10 



V. Khachatryan et al. [CMS Collaboration], JHEP 
1002, 041 (2010). 

V. Khachatryan et al. [CMS Collaboration], Phys. Rev. 
Lett. 105, 022002 (2010). 

K. Aamodt et al. [ALICE Collaboration], Eur. Phys. J. 
C 65, 111 (2010). 

K. Aamodt et al. [ALICE Collaboration], Phys. Lett. B 
693, 53 (2010). 

K. Aamodt et al. [ALICE Collaboration], Eur. Phys. J. 
C 68, 345 (2010). 

K. Aamodt et al. [ALICE Collaboration], Eur. Phys. J. 
C 68, 89 (2010). 

G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 
688, 21 (2010). 

G. Aad et al. [ ATLAS Collaboration], CERN ATLA S- 
CONF-2010-024, |http: / /cdsweb.cernxh /collection / 1 
ATLAS Conference Notes?ln=fr 

G. Aad et al.[ ATLAS Collaboration], CERN ATLAS- 
CONF-2010-031. 

G. Aad et al. [ ATLAS Collaboration], CERN ATLAS- 
CONF-2010-046. 

G. Aad et al. [ ATLAS Collaboration], CERN ATLAS- 
CONF-2010-047. 

G. Aad et al. [ ATLAS Collaboration], larXiv:1012. 51041 
[hep-ex] (submitted to New J. of Phys.). 
B. I. Abelev et al. [STAR Collaboration], Phys. Rev. C 
79, 034909 (2009). 

T. Aaltonen et al, CDF Collaboration, Phys. Rev. D 
79, 112005 (2009). 

T. Alexopoulos et at, E735 Collabora- 
tion, Phys. Lett. B 336, 599 (1994). 

C. Albajar et al, UA1 Collaboration, 
Nucl. Phys. B 335, 261 (1990). 

F. Abe et al., CDF Collaboration, Phys. Rev. D 41, 
2330 (1990). 

R. E. Ansorge et al, UA5 Collaboration, Z. Phys. C 
43, 357 (1989). 

R. E. Ansorge et al 
37, 191 (1988). 

F. Abe et al, CDF Collaboration 
1819 (1988). 

G. J. Alner et al. Phys. Rep. 154, 247 (1987). 
R. E. Ansorge et al, UA5 Collaboration, Z. Phys 

33, 175 (1986). 

W. Thome et al, Nucl. Phys. B 129, 365 (1977). 

N. Armesto et al, J. Phys. G 35, 054001 (2008). 

N. Armesto, in Quark Gluon Plasma 4, Eds R. C. Hwa 
and X. -N. Wang (World Scientific, Singapore, 2010), 
p375 l|arXiv:0903.1330[ hep-ph] ) . 

M. Mitrovski, T. Schuster, G. Graf, H. Petersen and 

M. Bleicher, Phys. Rev. C 79, 044901 (2009). 

E. Levin and A. H. Rezaeian, Phys. Rev. D 82, 014022 

(2010). 

E. Levin and A. H. Rezaeian, Phys. Rev. D 82, 054003 
(2010). 

Larry McLerran and Michal Praszalowicz, Acta Phys. 
Polon. B 41, 1917 (2010). 

C. Merino, C. Pajares, M. M. Ryzhinskiy, and 
Yu. M. Shabelski. larXiv: 1007.32"06l hep-ph] . 
P. Z. Skands, Phys. Rev. D 82, 074018 (2010). 
Andy Buckley, Hendrik Hoeth, Holger Schulz and 



UA5 Collaboration, Z. Phys. C 



Phys. Rev. Lett. 61, 



C 



Jan Eike von Seggern, PoS ACAT08, 112 (2009). 
[33] A. Moraes, C. Buttar and I. Dawson, Eur. Phys. J. C 

50, 435 (2007). 
[34] J. F. Grosse-Oetringhaus and K. Reygers, J. Phys. G 

37, 083001 (2010). 
[35] I. Kraus, J. Cleymans, H. Oeschler and K. Redlich, 

Phys. Rev. C 79, 014901 (2009). 
[36] A. B. Kaidalov and M. G. Poghosyan, Eur. Phys. J. C 

67, 397 (2010). 
[37] E. K. G. Sarkisyan and A. S. Sakharov, Eur. Phys. J. C 

70, 533 (2010). 
[38] Andreas Warburton for CDF and DO Collaborations, 

PoS HCP2009, 013 (2009). 
[39] M. H. Seymour, to appear in the Proceedings of Physics 

at LHC 2 010 (PHLC2010) , 7-12 June 2010, Hamburg 

Germany, larXiv:1008.2927l [hep-ph] . 
[40] R. Sassot, P. Zurita and M. Stratmann, Phys. Rev. D 

82, 074011 (2010). 
[41] Wei-Tian Deng, Xin-Nian Wang and Rong Xu, 

I arXiv: 1008.18411 [hep-ph] . 
[42] X. N. Wang and M. Gyulassy, Comput. Phys. 83 (1994) 

307; X. N. Wang, Phys. Rep. 280, 287 (1997). 
[43] S. E. Vance and M. Gyulassy, Phys. Rev. Lett. 83, 1735 

(1999). 

[44] B. Andersson, G. Gustafson, B. Nilsson-Almqvist, Nucl. 

Phys. B281, 289 (1987). 
[45] A. Capella, U. Sukhatme, C. -I. Tan et at, Phys. Rept. 

236, 225-329 (1994). 
[46] T. S. Biro, H. B. Nielsen and J. Knoll, Nucl. Phys. 

B245, 449 (1984). 
[47] M. Gyulassy and A. Iwazaki, Phys. Lett. B. 165, 157 

(1985). 

[48] F. Gelis, T. Lappi and L. McLerran, Nucl. Phys. A828, 
149 (2009); T. Lappi and L. McLerran, 
Nucl. Phys. A772, 200 (2006). 

[49] L. McLerran, J. Phys. G 35, 104001 (2008). 

[50] T. Sjostrand, Comput. Phys. Commun. 82, 74-90 
(1994). 

[51] X. N. Wang and M. Gyulassy, Phys. Rev. D 44, 3501 

(1991) . 

[52] X. N. Wang and M. Gyulassy, Phys. Rev. D 45, 844 

(1992) . 

[53] P. Levai and B. Muller, Phys. Rev. Lett. 67, 1519 
(1991). 

[54] T. Pierog, S. Porteboeuf, I. Karpenko and K. Werner, 
to appear in the Proceedings of 45th Rencontre de 
Moriond, QCD and High Energy Interactions, 13- 
20 March 201 0, La Thuile, Aosta Valley, Italy, 
larXiv:1005.4 526 [hep- ph], 

[55] A. Iwazaki, larXiv:0904.1449l [hep-ph], 

[56] J. S. Schwinger, Phys. Rev. 82, 664 (1951). 

[57] D. Kharzeev, E. Levin and K. Tuchin, Phys. Rev. C 
75, 044903 (2007). 

[58] N. Tanji, Annals Phys. 324, 1691 (2009). 

[59] R. Ruffini, G. Vereshchagin and S. S. Xue, Phys. Rept. 
487, 1 (2010). 

[60] N. Tanji, Annals Phys. 325, 2018 (2010). 

[61] T. D. Cohen and D. A. McGady, Phys. Rev. D 78, 
036008 (2008). 

[62] F. Hebenstreit, R. Alkofer and H. Gies, Phys. Rev. D 
78, 061701 (2008). 



11 



V. Topor Pop, M. Gyulassy, J. Barrette, C. Gale, 
Phys. Rev. C 72, 054901 (2005). 

V. Topor Pop, M. Gyulassy, J. Barrette, C. Gale, 
S. Jeon and R. Bellwied, Phys. Rev. C 75, 014904 (2007) 
(and references therein). 

N. Armesto et ai, J. Phys. G 35, 054001 
(2008); V. Topor Pop et al, ibid., pl5; V. Topor Pop 
et al., ibid., p 57. 

V. Topor Pop, J. Barrette and M. Gyulassy, 
Phys. Rev. Lett. 102, 232302 (2009). 

V. Topor Pop, M. Gyulassy, J. Barrette, C. Gale, 
X. N. Wang, N. Xu, 
Phys. Rev. C 70, 064906 (2004). 

G. Ripka, Lecture Notes in Physics, 639, 138, 2004 (Ed. 
Springer, Berlin, Germany). 

M. Cristoforetti, P. Faccioli, G. Ripka, and 
M. Traini, Phys. Rev. D 71, 114010 (2005). 
K. Nakamura et al. [Particle Data Group], J. Phys. G 
37, 075021 (2010). 

N. S. Amelin, N. Armesto, C. Pajares, and 

D. Sousa, Eur. Phys. J. C 22, 149 (2001). 
L. V. Gribov, E. M. Levin and M. G. Ryskin, Phys. 
Rept. 100, 1 (1983). 

G. Altarelli and G. Parisi, Nucl. Phys. B126, 298 
(1977). 

X. N. Wang, Phys. Rev. D 43, 104 (1991). 

E. Eichten, I. Hinchliffe, K. D. Lane and C. Quigg, Rev. 
Mod. Phys. 56, 579 (1984) [Addendum- ibid. 58, 1065 
(1986)]. 

D. W. Duke and J. F. Owens, Phys. Rev. D 30, 49 
(1984). 

D. Acosta et al., CDF Collaboration, Phys. Rev. D 65, 
072005 (2002). 

M. J. Tannenbaum, talk given at the Worshop on 
RHIC Par adigms, April 14-17,2010, Austin, Texas, 
larXiv: 1008. 15361 [nucl-ex] . 
B. Hippolyte, Eur. Phys. J. C 49, 121 (2007). 
B. Hippolyte, Eur. Phys. J. C 62, 237 (2009). 

H. Ricaud, A. Kalweit and A. Maire, J. Phys. G 37, 
094049 (2010). 

[82] R. J. Fries, B. Muller, C. Nonaka and S. A. Bass, Phys. 
Rev. C 68, 044902 (2003). 



[83 
[84 
[85 
[86 
[87 



V. Greco, C. M. Ko and P. Levai, Phys. Rev. C 68, 
034904 (2003). 

R. C. Hwa and C. B. Yang, Phys. Rev. C 67, 034902 
(2003). 

B. I. Abelev et al. [STAR Collaboration], Phys. Rev. C 
75, 064901 (2007). 

D. E. Acosta et al. [CDF Collaboration], Phys. Rev. D 
72, 052001 (2005). 

CMS Collaboration, CERN CMS Physics 
Analysis Summary QCD-10- 007 (2010), 
|http://cdsweb.cern.ch/record/1279344| Keith A. 
Uhner for the CMS Collaboration. larXiv:1012. 16051 
[hep-ex] . 

T. Sjostrand, S. Mrenna and P. Z. Skands, Comput. 
Phys. Commun. 178, 852 (2008). 
[89] T. Sjostrand, S. Mrenna and P. Z. Skands, JHEP 0605, 
026 (2006). 

[90] K. Werner, T. Hirano, I. Karpenko, T. Pierog, S. Por- 
teboeuf, M. Bleicher and S. Haussler, Nucl. Phys. Proc. 
Suppl. 196, 36 (2009). 

[91] K. Werner, I. Karpenko and T. Pierog, J. Phys. Conf. 
Ser. 230, 012026 (2010). 

[92] T. Alexopoulos et. al, E735 Collaboration, Phys. Rev. 
D 48, 984 (1993). 

[93] M. Banner et. al, UA2 Collaboration, Phys. Lett. B 
122, 322 (1983). 

[94] P. Antonioli for the ALICE Collaboration, to ap- 
pear in the Proceedings of Hadron Collider Sympo- 
sium (HCP2010 ), August 23-27, 2010, Toronto, Canada, 
larXiv:1010.3735l [hep-ex] . 

[95] V. V. Skokov and P. Levai, Phys. Rev. D 78, 054004 
(2008). 

[96] Peter Levai and Vladimir Skokov, PoS EPS-HEP2009, 
456 (2009). 

[97] A. Giovannini and R. Ugoccioni, Eur. Phys. J. C 36, 
309 (2004). 

[98] D. Kharzeev, Phys. Lett. B 378, 238 (1996). 
[99] A. K. Aamodt et al. [ALICE Collaboration], Phys. Rev. 
Lett. 105, 072002 (2010). 
[100] B. Z. Kopeliovich and B. Povh, Phys. Lett. B 446, 321 
(1999). 



